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We derive the spin Hamiltonian for the quantum spin liquid Na4lr30g, and then estimate the 
direct and superexchange contributions between near neighbor iridium ions using a tight binding 
parametrization of the electronic structure. We find a magnitude of the exchange interaction com- 
parable to experiment for a reasonable value of the on-site Coulomb repulsion. For one of the two 
tight binding parametrizations we have studied, the direct exchange term, which is isotropic, dom- 
inates the total exchange. This provides support for those theories proposed to describe this novel 
quantum spin liquid that assume an isotropic Heisenberg model. 

PACS numbers: 75.10.Dg, 75.30.Et, 75.50.Ee 



I. INTRODUCTION 

Na^^Os lias captured much attention since its 
discovery^ This insulator exhibits a large Curie Weiss 
temperature (650 K), yet does not magnetically order 
down to the lowest measured temperature. The reason 
for this is thought to be due to the strongly frustrated 
nature of the iridium lattice, which forms a hyperkagome 
network composed of corner sharing triangles (the hyper- 
kagome lattice being formed by replacing one of the four 
Ir sites in a pyrochlore lattice by a Na ion) . Most models 
that have been proposed to describe this material assume 
an isotropic Heisenberg model with an effective spin of 
1/2. This is somewhat of a surprise, given the distorted 
nature of the lattice and the strong spin-orbit coupling 
of the iridium ions. These issues have been discussed in 
depth in Ref. y. 

In a previous paper ^ we have calculated the ex- 
change constants based on a particular tight binding 
parametrization of the electronic structure, and found 
that the resulting spin Hamiltonian should be highly 
anisotropic. In the present paper, we revisit this issue 
by considering a more general tight binding parametriza- 
tion. We now include the residual crystal field splittings 
of the Ir 5d orbitals due to the octahedral distortions, and 
find that this corrects a major deficiency of the previous 
fit, which was an anomalously large value for the t s dd hop- 
ping. As a result, the more general parametrization leads 
to the exchange interaction being dominated by the di- 
rect exchange between Ir ions, and as a consequence, we 
find an approximately isotropic Heisenberg model. We 
also find the observed magnitude of the exchange for a 
reasonable value of the Coulomb repulsion. 

In Section II, we provide a microscopic derivation of 
the exchange Hamiltonian, finding a few differences from 
previously published results. In Section III, we describe 
the tight binding parametrizations of the electronic struc- 
ture, and then in Section IV, the resulting exchange con- 
stants as a function of the Coulomb repulsion. In Section 



V, we summarize our findings. 



II. EXCHANGE COUPLINGS 

In Na4lr 3 08, the t2 g 5d manifold contains one hole per 
iridium site (the e g levels are empty), and the oxygen 
2p levels are filled 4^ This single hole sits in a half-filled 
doublet level due to spin-orbit coupling, thus motivating 
an effective S=l/2 exchange model. This picture, which 
is supported by our electronic structure calculations, can 
be exploited to calculate the exchange couplings following 
Ref. @. 



A. Microscopic Hamiltonian 

In our own derivation, we exploit the work of Refs. 
which calculated the exchange couplings for cuprates in 
the presence of spin-orbit coupling, following the earlier 
work of Koshibae et alJ- referred to in Ref. 3. We want 
to derive the exchange couplings generated to fourth or- 
der in the hopping of Ir holes in the ti g complex. Hop- 
ping can be direct between the Ir 4+ ions or via the 



p-orbitals of the O ions. We denote by 



d 

e l,2,3 



the 



Kramers-degenerate energy levels of the ti g complex with 
ef > ef , e%. These splittings are mostly due to the spin- 
orbit coupling, which acts to form a lower quartet and an 
upper doublet, but there is a contribution as well from 
residual crystal field splittings resulting from the low site 
symmetry of the iridium ion (noting that the space group 
is cubic).— >^ The vacuum state is denoted as |fi) where all 
three levels are fully occupied. The ground state is then 
generated by linear combinations of the form 



|0o) = span{ct 3t |0),ct 34 jn)} 



(1) 



where f>4- characterizes the Kramers degenerate states, 
i is the site index, and is the creation operator for a 
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hole in the t^g complex. According to the Goodenough- 
Kanamori rules, the strongest contribution to the ex- 
change coupling results from two half-occupied orbitals, 
so we may focus only on the hopping between the eg or- 
bitals. The Hamiltonian acting on the ground state is 
then 



H = H + Hhop 
with on-site Hamiltonian 

TT _ \ " d j I d \ ~* t t 



(2) 
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Up j i 
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(3) 



where c' (c) and (p) are creation (annihilation) oper- 
ators for holes. Indices i and k describe sites of the Ir 
ion and its nearest-neighbor O ions, m, n and a, b denote 
the Ir tig and O p levels, and a, a' are spin indices. The 
relevant hopping between nearest neighbor Ir and O sites 
is 

-Hhop = (^33)^' C i3o- C j3cr' 

+ E ((*S)-'pL^+ii.c.) (4) 

zfcatTcr' 

where the spin dependent matrix elements for hopping 
between Ir, and Ir and O ions, are 



1>33,W — ^33"cr(7' T ^33 "crcr' 
(^a3) CTCr ' — *(j3<W' + C^ 3 • CT^/ 



(5) 
(6) 



Note that we use a formulation in terms of holes. But to 
compare with the results of Ref. 3j, we will insert for the 
energies of virtual states that of the electrons and not 
the holes. 



B. Exchange couplings from perturbation theory 

The Kramers degenerate ground state is split due to 
virtual hopping processes. An effective spin Hamiltonian 
Hs can be derived from perturbation theory in H\ lop . Ac- 
counting only for interactions between neighboring spins 
it is of the general form 



(i0) 



(7) 



Here (ij) indicates a sum over nearest neighbor sites and 
Mid =J2 J P q(hj)S p {i)S q {j) (8) 

pq 

with spin-1/2 operators S (p, q — x, y, z). In what follows 
we want to derive expressions for J pq {i,j) up to fourth 



order in the hopping. For this it is convenient to intro- 
duce 





00' 


(9) 


rppd 


= E (^3 W'PL CT c i3<x' 
acrcj' 


(10) 




acrcr' 


(11) 


H hop = E 


^f+E(nf + (^)! fe ) 

ki 


(12) 



so that 

The energy of the ground state |0o) is se t to zero. 

1. Direct exchange 

The lowest order contribution results from direct hop- 
ping between Ir ions. There are two contributions result- 
ing from the back and forth hopping of holes on Ir sites 
i and j, respectively. Both processes give identical con- 
tributions, so we can restrict ourselves to the hopping of 
a hole at site i and multiply its contribution by a factor 
of two, 



U?)=-2{<h\Tgj£r$t\<h) 



dd 1 rpdd\ 



(13) 



As mentioned above, we want to compare with the find- 
ing s of Ref. i and therefore we give the energy of the in- 
termediate state with two holes on Ir site j in terms of the 
electron's energy. The latter is — 4Ud + 5Ud = Ud and re- 
sults from the difference in reduced (increased) Coulomb 
interaction on the site where a hole is inserted (removed). 
Inserting Tff from Eq. 9 we get 



(2) 



jj E/ (^33)o"40- 3 (£33)0-2^1 



(Tl <T2 <T3 <T4 



{Mcl-ia^j^A^ ^ I'M (I 4 ) 

We now employ the identity ( for notational convenience 
we suppress the unit matrix in ^1) 



C i3cri C i3o 

to rewrite this as 



'h + Si ■ a) 



(15) 



nfl = 



, : l ^tr(4(I + Sr <r)4(5 + V *)) (16) 

where the trace is over spin-indices and we neglected a 
quadratic term cc 1 which does not contribute to the ef- 
fective spin Hamiltonian in Eq. 8. Terms involving the 
factor 1/2 lead to contributions which also do not con- 
tribute to Eq. 8. Neglecting these we find 

«g ^r(r-AS,-n:,!: :i iS,-,:) (17) 

This expression is a factor of two larger than that quoted 
in Ref. U 
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2. Superexchange via oxygen 

We now turn to processes involving intermediate oxy- 
gen and states arising from fourth order hopping. Let us 
start with two observations. First, only those processes 
where the hole hops from one Ir via the O to a differ- 
ent Ir give spin dependent contributions, i.e. processes 
where the hole hops to O and then returns to the same 
Ir it started from do not contribute to Eq. 8 and will be 
neglected in the following. Second, following Ref. 0, we 
notice that there are two qualitatively different processes 
in which fourth order hopping contributes to Eq. 8. In 
the 'consecutive channel' (cc) a hole hops from an Ir ion 
to the O, then to the second Ir, and finally back to the 
first Ir via the same or a different O. In this channel there 
are two holes on an Ir ion in the intermediate state. The 
consecutive channel cc gives an identical contribution as 
the direct term, however with an effective hopping am- 
plitude £33 — > T33. In the second 'simultaneous channel' 
(sc) a hole hops from one of the Ir ions to the O, and 
then a second hole hops from the second Ir to the same 
or a different O. Afterwards the holes return to the Ir, 
i.e. back to the ground state in which one hole occupies 
each Ir. In this second channel, two holes simultaneously 
occupy O as intermediate states. 

a. Consecutive channel cc: Let us first turn to con- 
tributions from the consecutive channel (cc) and look at 
a hole at Ir site i. Hopping of the hole from site i to j via 
an O ion at site k results in a state with two holes 
occupying ti g levels of Ir at site j and no hole on Ir site 
i, 

H \<fi%) ^%-Lt£Vo> (18) 

The (electron) energy of the intermediate state after the 
first hop is a sum of four contributions. Removing the 
hole from the ti g level and inserting the hole in the a- 
level of the O ion at site k gives ef — e p ka for the on-site 
energies and 5(U~d — U p ) from the increased and reduced 
Coulomb interaction between the electrons on Ir and O 
ions, respectively. That is, the energy of the intermediate 
state after the first hop is 

4t = 4-4a+HU d -U p ) (19) 

Inserting Eq. 10 then gives 

Ho\<t>%)= E i(€w 2 (&) nn 

X c ]3<T3-P^2-PL<t 2 c »3<tiI0o) 
= E -Jd^la)^ ©^.r, c] ia3 C aai \(t>o) 

(20) 

Introducing the effective hopping between Ir sites 

4t=J2^f (21) 

ak ^ka 



which is a matrix in spin space (summation over inter- 
mediate spins is implicit), this is rewritten as 

H °\<t%) = E (^W^ks^ \<h) (22) 

Eq. 22 corresponds to the intermediate state T^ d |(/>o) in 
the direct exchange hole pathway Eq. 13 with the hop- 
ping matrix elements T33. Therefore division by the (elec- 
tron) energy Ud of the intermediate state Eq. 22 and ap- 
plication of the second hop which returns the hole to Ir 
site i, we find 

^f ) =|-tr(4(S J . ( r)^(S t .a)) (23) 

where the factor of two accounts again for the fact that 
an identical contribution results from a process in which 
the hole starts at site j. 

b. Simultaneous channel sc: In the 'simultaneous 
channel' sc both Ir-holes from sites i and j hop onto O 
ions in the intermediate state |</>| c ). This can occur in a 
process in which the first hole at site i and then the hole 
at site j hops, or in the reversed order. The first hop in 
this process results in the state 

k 

= E/ i^az)cr2<yiVkaa 2 C i^l 
ak(T\ <T2 

+ (^3)^1^3^)100) (24) 

The energy of the intermediate state \4>\ c ) is again e]L in 
Eq. 19. Adding then the second hole on the same or a 
different O ion, we obtain the intermediate state H^cf)^), 
i.e. 

WS C > = E ( T ujf T v + T «w T ^ ) l0o) (25) 

An explicit expression for this intermediate state is given 
in Appendix A, and we here only remark that its (elec- 
tron) energy depends on whether both holes are on the 
same O ion or not. It has the non-interacting contribu- 
tion 2eg — e p ka — ef b , an interaction contribution 5Ud + 5Ud 
from removing holes at Ir sites i and j, and the interac- 
tion contribution —9U p if both holes are on the same O 
ion and — 10U P in case they are not. This can be sum- 
marized as 

^ d + 4t + U P 5 kl (26) 

The return to the ground state can occur again in two 
ways, i.e., a hole returns first to Ir at site i and then 
the second returns to site j, or in the reversed order. 
The overlap of the resulting state with the ground state 
defines the effective Hamiltonian 
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(27) 



The calculation of Eq. 27 is identical to the one for 
cupratcs in Rcfs. 5,6. For completeness we give details in 
Appendix B, and here merely state that neglecting spin 
independent contributions, the effective Hamiltonian of 
interest is 



(4sc) 



E 



balk c lb 



pd . pd 



pd pd 
> e /6 e ka, 



UpSkl 

xtr (tgtgiS^^ttt^iSi-a) 



(28) 



Finally, summing the superexchange contributions from 
both channels the superexchange contribution is 



balk 

where we introduced 



(29) 



(2) 

to find that T-L\ ■ is of the form Eq. 31 with 



7(2) _ pi 



D (2) 

r 



► (2) 
ij 



4i 
4 Z^- 7 



(7^ V 



33^33 



33^33 
ij -tji 



) 



( C33 C 33 + C33C33 C33 • C 33 1 



(34) 
(35) 
(36) 



To do the same calculation for the superexchange Eq. 29, 
let us first rewrite products of hopping matrices 



l Jj _ (Z: 



L 3b i b3 



(t 



36 



63 



^63 



— i il flj j_ c il r 1 ^ 

— c 3b c 63 ^36 ' *^b3 



^(^36 X ^63) + ^36^ 



b3 



where in the last line we introduced 



(37) 



Jk 



1 



JP<1 _i_ ^P<i 
fc /6 



::+UpS kl 



pd 
6 ka, 



(30) 



Note that the expression for g differs from that of Ref. [3j. 
This difference was previously noted in Refs.06,8 in con- 
nection with the similar expression in Ref. [jj. 



C. Exchange Constants 

We next want to bring Eq. 17 for the direct exchange 
into the form 



Hi 



JSi ■ Sj + D lJ ■ (Sj x S,0 + Si ■ T -8j (31) 



To this end we insert Eqs. 5 and 6 into Eq. 17 



■1/(2) _ 
^,3 - 



(32) 



fil pj 



3b°b3 ^36 ' "^63 

-•//-lii v /~i!j \ 1 IS! /-i!j , /-<i! rZj 
l[\> 3b X ^ &3 J + S 3b ^ 63 + ^3b l b3 



J 3b ~ ^63 

We can now apply again Eq. 33 to find 



(38) 
(39) 



l/W _ 
H *,3 ~ 



E^tr((4+vg-a) (S, -a) 
x ( S *?+v£.<r)(S,-<7) 



=JS I • S 3 +D' iJ • (S, x Sj) + S, • T -Sj (40) 



with 

J(4) 



= 2E 



Jb Ik ka 
s ij!Jba ji 



balk 



D 



<-, • /",J6 Ik ka 

-'1^ { v ij9ba s ji 



balk 



r ^ = 2 E r^^. 



bulk 



\ku 



lb Ik ka 
s ij9ba v ji 



4— lb — ^.fca 
ji9ba * 



(41) 
(42) 



y tj.Jba v ji 



(43) 



We then employ the identity (checked in Appendix B) 



and <7{,^ given in Eq. 30. 



~tr ((A! + Bx • <r)(S 4 • a)(A 2 + B 2 • a)(S 3 ■ a)) 
= A 1 A 2 S i -Sj+i (A 2 Bi - AiB 2 ) • (S, x Sj) 

+ Si- (BiBa+BaBi-Bi-Bal) -Sj (33) 



III. TIGHT BINDING PARAMETRIZATION 

In Ref. 0, we performed local density approximation 
calculations for the electronic structure of Na^^Os, and 
then performed a tight binding fit within a Slater-Koster 
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formalism. To understand the tight binding fit, we need 
to step back and take a look at the electronic structure 
of this material.—^ Na^^Os is composed of IrOg octahe- 
dra, and as with many transition metal oxides, these oc- 
tahedra are distorted, with just a C2 symmetry axis pre- 
served. Moreover, there are two types of oxygens, with 
four of the six around an Ir ion being of one type (02), 
the other two of the other type (01). Note that the two 
01 ions are not related by an inversion like in cuprates, 
but are related by a ir rotation about the C2 axis (see 
Fig. 1). A minimal tight binding model restricted to 
near neighbors only would then consist of hoppings of Ir 
to 01, Ir to 02, 01 to 01, 02 to 02, 01 to 02, and Ir 
to Ir. Because of the distorted nature of the lattice, not 
all distances between given atom types are the same. To 
reduce the number of fit parameters, we then assumed a 
typical inverse fourth power dependence of the hopping 
integrals with distance for a given atom combination. In 
practice, this affects only the O to O hoppings, and we 
note that an inverse fourth power behavior was indeed 
found for Ir02<^ The net result is that we need 18 tight 
binding parameters: energies of t2 g (et 2s ) ancl e g ( e e g ) 
5d orbitals on Ir, energies of 2p orbitals on 01 and 02 
(eoi, £02) d-d hoppings (t% d , q d , t s dd ), d-p hoppings (t% p , 
t^ p for Ir-01 and Ir-02), p-p hoppings (i£ p , i£ p for 01- 
01, 02-02, 01-02), and the coefficient of the spin orbit 
splitting for the Ir 5d orbitals, A (i.e., A 1 • s). 

These parameters are then used to evaluate the var- 
ious elements of the secular matrix that generates the 
eigenvalues. 10 The diagonal elements are simply given 
by the various on-site energies, and the spin-orbit cou- 
pling matrix elements in a crystal field basis com- 
patible with that used in Ref. [13 can be found in 
Ref. [U. The off-diagonal matrix elements are of the form 
tai,bj(l> m i n)e lk '( rj ~ r ^ where t is the hopping integral be- 
tween orbital a on site i and orbital b on site j, and I, m, n 
are the three direction cosines between the two sites at 
rj and Tj. In our case, the resulting secular matrix has 
dimension 312 (there are four formula units in the unit 
cell, i.e., there are 120 Ir 5d orbitals and 192 O 2p orbitals 
in the unit cell when spin-orbit is included) . In consider- 
ing off-diagonal elements of the secular matrix, note that 
each iridium ion is surrounded by six oxygens and four 
other iridium ions (Fig. 1), and each oxygen is surrounded 
by twelve other oxygens. Once the secular matrix is set 
up, then the various tight binding parameters are iter- 
atively adjusted to achieve an optimal fit to the band 
structure eigenvalues. The function being minimized is a 
sum of the squares of the differences of the tight binding 
eigenvalues from the first principles ones of the electronic 
structure calculation. The minimization was performed 
using Powell's method, 12 We first fit the calculation with- 
out spin-orbit using 42 eigenvalues (the bottom and tops 
of the 01, 02, and e g complexes, as well as all 36 bands 
of the t2 g complex) at each of the four symmetry points 
of the simple cubic Brillouin zone. As an initial start to 
the minimization, we used previously derived tight bind- 
ing parameters for Ir02, 9 ' scaled (by an assumed inverse 




FIG. 1: (Color online) Superexchange pathway between two 
Ir ions as marked by the arrows. The Ir ions are the large 
(red) spheres and the oxygen ions the small ones. The lower 
left (green) pathway is via an 02 (blue) ion, the upper right 
(purple) pathway via an 01 (brown) ion. 

fourth power dependence on distance) to the Na^^Os 
lattice. The starting values for the on-site energies were 
estimated from the orbital decomposed density of states 
of the band calculation. After minimization, the result- 
ing tight binding parameters were then used as input to 
fitting the calculation with spin-orbit, now involving 78 
eigenvalues (the bottom and tops of the 01, 02, and e g 
complexes, as well as all 72 bands of the t2 g complex), 
again at the four symmetry points. For the initial start, 
the spin-orbit coupling parameter, A, was estimated from 
the splitting of the j = 3/2 and j = 5/2 eigenvalues of 
the band calculation. After minimization, the resulting 
fit gave a good reproduction of the energy bands, along 
with the Fermi surface (Figs. 8 and 9 of Ref. 0). This is 
non-trivial, given the low site symmetry of this lattice, 
the size of the secular matrix, and the large 18 parameter 
function space. 

This fit, though, does not take into account the residual 
crystal field splitting of the 5d orbitals due to the distor- 
tions of the octahedra. These splittings, although some- 
what obscured by hybridization, appear to be present 
(Fig. 3 of Ref. [H), and could potentially be of importance. 
Including them increases the number of tight binding pa- 
rameters to 21 (with two more parameters needed for 
the tig manifold and one for the e g one). A complete 
description of these splittings are complicated because 
they involve a number of crystal field potential terms, 
V20, V21, V22, V i0 , V41, V i2l V i3l V44, whose coefficients are 
difficult to estimate from first principles. Rather, we used 
a simplified approach where we ignore the small coupling 
between the t2 g and e g orbitals. For a C2 axis along 
(1,1,0), the t2 g diagonal elements of the secular matrix 
would be (relative to et 2 ) of the form -2ci for xy, and 
+ci for xz and yz, with an off-diagonal matrix element 



TABLE I: Tight binding hopping parameters in eV from the 
21 parameter fit. The on-site energies are eoi = -6.4241, 602 
= -3.9141, e t2g = -1.7230, e eg = 0.6619, with the spin-orbit 
coupling A = 0.5797. The residual crystal field splittings of 
the cubic levels on the Ir sites are denoted as a [i = 1, 2, 3). 
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Ci 


0.0324 


-0.3839 


0.2965 



±C3 between xz and yz, the sign depending on whether 
the C2 axis is along (1,1,0) or (1,-1,0). Diagonalization 
of this sub-matrix leads to two even symmetry states and 
one odd symmetry state relative to the C2 axis. We note 
that this simplified form ignores the smaller off-diagonal 
matrix element between xy and xz, yz which would cou- 
ple the two even symmetry states. For the e g sub-matrix, 
the diagonal matrix elements (relative to e eg ) are -2c2 for 
x 2 — y 2 and +2c2 for 3z 2 — r 2 . Diagonalization of this 
sub-matrix leads to one even symmetry and one odd sym- 
metry state. As the C2 axis is rotated from one iridium 
site to the next, these t^g and e g sub-matrices in turn 
must be rotated (leading to off-diagonal terms between 
x 2 — y 2 and Zz 2 — r 2 ). 

The resulting 21 parameter tight binding fit had about 
a 25% smaller RMS error than the 18 parameter one (the 
band dispersion from the fit is plotted in Fig. 2). On the 
other hand, the solution space is more complex than in 
the 18 parameter case, and as a consequence, it is difficult 
to estimate how well the minimization routine has suc- 
ceeded in finding an optimal solution. In particular, the 
Ci parameters (Table I) lead to effective level splittings 
for the t2 g and e g orbitals that do not seem to correspond 
well to those indicated by the band structure calculation 
(Fig. 3 of Ref. 0). Moreover, the Fermi surface is some- 
what degraded relative to the one of the 18 parameter fit 
(Fig. 9 of Ref. 0). On the other hand, the 21 parameter 
fit corrects a major deficiency of the 18 parameter one, in 
that t s dd , which was anomalously large in the 18 param- 
eter fit (0.1545 eV), is now far more reasonable (0.0521 
cV). This, and other differences in the tight binding pa- 
rameters (comparing Table I with Table II of Ref. |j), 
turn out to have a qualitative impact on the exchange 
constants, as we will see in the next Section. 

IV. EXCHANGE CONSTANTS 

The exchange constants are derived by taking the tight 
binding parameters discussed in the previous Section and 
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FIG. 2: Energy bands (t2 g spin-orbit 'doublet') near the Fermi 
energy (Ef) from the 21 parameter tight binding fit to the 
electronic structure of Na^^Os- The horizontal line marks 
E F . 



inserting them into the expressions derived in Section II 
(with = 0) . We remind that the functional form for the 
spin Hamiltonian is listed in Eq. 31. As stated previously 
111 Ref. 1, the direct exchange term only contributes to 
the isotropic term J in Eq. 31; the other terms vanish, 
even for the distorted lattice. This is a consequence of 
the fact that the spin-orbit coupling leads to a 'doublet' 
of states around the Fermi energy (the 24 bands shown 
in Fig. 2) which to a good approximation are formed 
from linear combinations of xy, xz, and yz orbitals with 
equal weights. The resulting contribution to J is of the 
form 4*3/17 where t d = \t° dd + + & S M (this differs 
by a factor of two from Ref. 4, as already noted after 
Eq. 17). On the other hand, the superexchange terms 

•H- 

contribute to J, D and T. We denote the components 
of the diagonal contributions to the exchange by Ji = 
J + Di + Ta (i = x,y,z). These values are plotted in 
Fig. 3 as a function of Ud (Coulomb repulsion on the Ir 
sites) for both sets of tight binding parameters. For the 
purposes of these plots, U p (Coulomb repulsion on the O 
sites) was set to zero. 

For the 18 parameter fit, J x becomes equal to the ex- 
perimental value of 28 meV— for a value of Ud of about 
1.1 eV (Fig. 3a). In that context, although the value of 
Ud is not known for Na^^Os, we note that for the re- 
lated perovskite, Sr2lr04, the optical gap is 0.5 eV, and 
it is known from LDA+U simulations that a Ud of about 
2 eV is needed to reproduce this gap. 1 ^ Interestingly, J y 
and J z significantly differ from J x , indicating that the 
predicted spin Hamiltonian from this tight binding fit is 
strongly anisotropic, as we previously remarked. 4 

We can contrast this with the 21 parameter fit, shown 
in Fig. 3b. J x reaches the experimental value of 28 meV 
for a Ud of about 1.5 eV, which is close to the antici- 
pated value of 2 eV. More interestingly, J y and J z are 
close in value to J x . This indicates a far more isotropic 
spin Hamiltonian, which is in support of various theo- 
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FIG. 3: (Color online) Diagonal elements of the exchange in- 
teraction as a function of the on-site repulsion (Ud) on the 
iridium sites, assuming zero repulsion (Up) on the oxygen 
sites; (a) is from the 18 parameter tight binding fit and (b) 
from the 21 parameter fit. The horizontal dashed line is the 
experimental value for J. 1 



FIG. 4: (Color online) Diagonal elements of the exchange 
interaction as a function of the on-site repulsion (U p ) on the 
oxygen sites, with a repulsion (Ud) of 1.5 eV on the iridium 
sites; (a) is from the 18 parameter tight binding fit and (b) 
from the 21 parameter fit. The horizontal dashed line is the 
experimental value for J.— 



ries for this material. For completeness, we list all the 
coefficients of the spin Hamiltonian in Table II for the 
21 parameter fit with Ud = 1.5 eV. We note not only 
the effective isotropy of Ji, but also the much reduced 
value of the Dzyaloshinski-Moriya interaction compared 
to what was previously indicated in Ref. 0. This differ- 
ence is mainly due to the much smaller value of Ud (0.5 
eV) assumed in the previous work. On general grounds, 
we note the dominance of Jd in Table II compared to 
the other exchange terms, in particular, the large ratio 
of Jd to J s . This is in contrast to the well known case 
of cuprates, where the superexchange term is dominant. 
This difference can be attributed to the 90 degree Ir- 
O-Ir bond present in this material compared to the 180 
degree Cu-O-Cu bond found in the cuprates. Other qual- 
itative differences between the 90 and 180 degree cases 
have been emphasized in the recent work of Jackeli and 
Khaliullin,ii 

Finally, we show in Fig. 4 the dependence of the ex- 
change constants on U p . Its effect is to cause increased 
anisotropy. However, this is more pronounced for the 18 



parameter fit than the 21 parameter one. At larger values 
of Up than those shown in Fig. 4, there are divergences 
that are associated with zeros of the denominators enter- 
ing g (Eq. 30). This corresponds to intermediate states 
that are low in energy, meaning that perturbation theory 
is no longer valid. Obviously, for very large values of U p , 
the superexchange contributions disappear, leaving only 
the isotropic direct exchange term. 

From the above, it is obvious that the exchange con- 
stants are very sensitive to the tight binding parametriza- 
tion, and also the values of the various Coulomb repul- 
sions. We also note that higher order processes are ig- 
nored, for instance, sixth order superexchange processes 
involving additional hoppings between the two O ions 
connecting two Ir sites which are known to play a role 
for the honeycomb lattice found in Na2lr03. 15 Longer 
range hoppings on the iridium sublattice could also be of 
importance as wellJ^ 



TABLE II: Exchange constants in meV from the 21 parameter 
tight binding fit. The spin Hamiltonian is given in Eq. 31. 
Quoted are values where Ir site m is along an (0,1,-1) direction 
relative to Ir site n. D are the Dzyaloshinski-Moriya, and T 
the anisotropic superexchange terms (in the second row for 
r, the jk refer to the parenthesis, xy, etc.). In addition, the 
direct exchange (isotropic) is Jd = 24.9, and the isotropic 
superexchange term is J a — 2.2. The last row, Ji, is the total 
exchange for the diagonal components (Jd + J s + r«). The 
assumed value of Ud is 1.5 eV (with U p assumed to be zero). 



i (jk) x (xy) y (xz) z (yz) 

Di 3.4 0.4 -0.3 

V lt 1.3 -1.3 -1.3 

I> 0.3 -0.2 -0.0 

Ji 28.4 25.8 25.8 



V. SUMMARY 

Wc have found that an approximately isotropic Heisen- 
berg model can be motivated from a tight binding 
parametrization of the electronic structure of NatlrsOs, 
with the experimental value of J reproduced for a rea- 
sonable value of the Coulomb repulsion, Ud- To obtain 
this result, it was important to account for the resid- 
ual crystal field splittings of the Ir 5d orbitals due to 
the octahedral distortions. Our findings are obviously of 
some importance in regards to models for this quantum 
spin liquid, since anisotropy acts to stabilize long range 
magnetic order.— The large value of the exchange is of 
much interest, since a large J appears to be associated 
with the unusual properties of cuprates, including their 
d-wave superconductivity^ 
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Appendix A: Intermediate states in the superexchange pathway 

For self-containedness this Appendix provides intermediate steps of the superexchange calculations similar to those 
found in Rcf. H The explicit expression for the intermediate state in the simultaneous channel of the superexchange 
pathway, Eq. 25, is a sum of the two contributions 



wr> =E I ?iT r « = -EE E E ^(^^(^j^pLpU^^i^,) (ai) 

Ik ® ba Ik 0403 faeri £ ka 

tfoic 6 ) =J2 T H^-nf\<Po) - -EE E E ^(^^(^j^pUpU^^i^) (a 2 ) 

Ik ba Ik 0403 0201 ka 

where the minus sign results from exchanging operators p' and c. Upon relabeling indices a b, k •<-> I, u<z •<-» 174, 
and fii O- 03 we add up both contributions, resulting in the intermediate state 



wr> = -EE E E (4 + 4r) ($w£wLj»U c . 

fca ifc 0403 o 2 oi \ e i6 e fca / 



•j3o- 3 Cj3o-i 1 4>o) (A3) 

and the state |</>| c ) has the energy Eq. 26. The return to the ground state can occur in the two processes described 
in Eq. 27. In the intermediate state i?o| </>§£) = Em^lLl^f) 01 tne nrs t process one hole has returned to site i, 

zj \J.sc\ (^3c l )o 6 05 (^63)040-3 (^3)0-20-1 ( 1 lit t t 1 / \ / A /i \ 

ffo|03i) = - 2^ 2^ <L 2^ pd, P d \tt x + -Jd) c a<y 6 P™c* 5 P[ b<j y ka<T2 c^ 3 c aai |0o) (A4) 

cba mifc CT6O50-403 o 2 ffi £ Zb + £ fca ' U p°kl \e ;f) £ fea / 

Using fermion anti-commutation relations and keeping only those contributions that are not annihilated upon acting 
on the ground state, Pmccr & Pi br74 P ka<J2 = 3 m idbcSa 5 a 4 p{ a(72 ~ S km S ac 5 a5r72 pj ba . 4 , which inserted into Eq. A4 shows that 

i?o|^) = J ffo(|03r) + IC' 6 )) where 

tj \isc,a\ V 1 (^36)0604(^63)0403(^03)0201 / I 1 \ + t I / \ /lr\ 

W 3 * ) = 2^ 2^ 2^ 2^ pd , ^~771 1 ^2 + -pd J PLo2 C I 3 o 6 9303^301 |0o) (A5) 

ba Ik cr 6 a 4 a 3 o 2 oi € lb + £ fca + ^P ^^ 

04 03 



#oic 6 > = -EE E E 



1 




% 


pd)' 

"ka/ 


(± 


1 


\ e pd 

\ e lb 


+ -p~d 

fc fea 



States |</>3i' a ) and ) have the energies e p ka and e p lb , respectively. Acting with a second hopping Hamiltonian the 
second hole returns to its ground state 

E(rppd\\ \ r L sc - a \ _ (4a)o8Q2 (4kig (^o3) cr 4Q3 ftal)^^ ( 1 I \_ \ J J „ „ U \ 

K 1 IjmWzi 1—2.^2-^ 2-^1 / pd , pd ,7-7 r \ pd |»i + pd I C j3a s c i3o 6 c j3o 3 C l3(Tl |^ / 

m ba Zfc J 8 iT6fr4ii3»2iTl v C ii> ' € ka ' U P°kl ) e ka \ C lb C ka / 

(A7) 

O8O4 (^3a)o6Q2 (^03)0403 (^3)0201 I 1 1 1 f t U \ 

V 1 )jm\93i /— 1^,2-^ 2.^ 2^ , pd pd jj c , pd \ pd + pd I C i3o 8 ^306^303^3^ | <po) 

m ba Ik osanoioztj^ox \ e lb ' e ka ' U P° kl ) € lb \ e lb £ ka / 

(A8) 

Processes in which the holes return to the ground state in reversed order have the intermediate state \4>3j) in which 
Ir ion at site j and i are in the ground and excited states, respectively. This intermediate state is again a sum of 
two contributions which are taken to the ground state \4>o) upon acting with the second hopping Hamiltonian. The 
final state corresponds to Eqs. A7 and A8 upon exchanging i f> j in the first and second matrix elements (i.e., those 
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returning the holes to their ground states) and the cr's. Upon relabeling of spin indices these four contributions to 
the superexchange pathway from the simultaneous channel can be combined into a sum of two contributions, whose 
overlap with the ground state is 



m(EE E E 



fjks /.it s f Jj \ /,ki • 



pd . pd , jt c 
ba Ik <7 8 a- 6 (j 4 (7 3 a 2 <?l € lb ' € ka ' u p°kl 



^gjggffg (^3b)g6g4 ftb3) CT 4<T3 (^3)^1 / 1 1 \ J t „ „ \u \ ( A Q\ 

1^1^ 1^ 1^ pd , pd , rr , \^d + -pd ] C i3<T 8 C j3<7e C j3<T3 C i3<7 1 I |0 O ) (A9J 



e». -he 



1 


±\ 


-pd 

% 


e pd J 

t kaJ 


f 1 


1 


\ t lb 


+ -Jd 

t ka 



Permuting some of the c, and applying identity Eq. 15, this can be summarized in the more compact form 

< sc) = E E pd- f 4? + 4r) tr (i + Si « (i + * - 

ba ifc % + 3a+ U P*W \4b 4a J V 

+ E E gT pd* f + -i) (44 (I + Sj - a)) tr (I + S, - a)) (A10) 

ba ife % +4a + U P S kl Vlb 4a) ^ ' 

where the trace is again in spin-space. The second term as well as contributions from 1/2 in the first term lead to 
spin independent contributions. Neglecting these, we arrive at the effective Hamiltonian Eq. 28. 

Appendix B: Traces over Pauli matrices to find 

Starting out from the expression 

tr ((Ax + Bi • a)(Si ■ a){A 2 + B 2 • a)(Sj ■ a)) 

= S*S]A x A 2 tt (aV) + S*S T A lB [tr (aVa™) 1 {{Ai + Bi . ^ . ^ + ^ . ff)(s . . a)) 



2 

= AxAiSi -S 3 +i (A 2 Bi - AiB a ) • (S 4 x S 3 -) 



5?5f A 2 Bf tr (a k a l a m ) + Sf SJ 1 B? B l 2 tx (a k a l a m a r - 

(Bl) 

■ , + Si ■ ( BiB 2 + B 2 Bi - Bi • B 2 1 ) • S 7 - (B5) 

we use the identities 1 V r ^ y J ^ ; 

tr (dV) = 24; (B2) 
tr (ctVct™) = 2ie fcim (B3) 
tr (oV<7"V n ) = 2(S k i5 mn - 5 km 6 ln + 5 kn Si m ) (B4) 
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